Periodic alternating 0, 7r-junction structures as realization of (p- Josephson junctions 
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We consider the properties of a periodic structure consisting of small alternating 0- and n- Joseph- 
son junctions. We show that depending on the relation between the lengths of the individual junc- 
tions, this system can be either in the homogeneous or in the phase-modulated state. The modulated 
phase appears via a second order phase transition when the mismatch between the lengths of the 
individual junctions exceeds the critical value. The screening length diverges at the transition point. 
In the modulated state, the equilibrium phase difference in the structure can take any value from 
— tv to tv (^-junction). The current-phase relation in this structure has very unusual shape with 
two maxima. As a consequence, the field dependence of the critical current in a small structure is 
very different from the standard Fraunhofer dependence. The Josephson vortex in a long structure 
carries partial magnetic flux, which is determined by the equilibrium phase. 



The minimum energy of Josephson junction (JJ) usu- 
ally corresponds to the zero phase difference of the super- 
conducting order parameter^. However, long time ago it 
has been predicted that the J Js with ferromagnetic inter- 
layer (S-F-S junctions) may have ground state with phase 
difference equal to ir ("7r-junctions")2. Only recently this 
prediction has been experimentally verified^*^ and now 
the controllable fabrication of S-F-S 7r-junctions becomes 
possible. Another realization of 7r-junctions became pos- 
sible due to the d-wave symmetry of the order parameter 
in high-temperature superconductors&t Recent exper- 
iments on YBa2Cu30y-Nb ramp long JJs of zigzag ge- 
ometry demonstrated that such junctions are composed 
of alternating facets of 0- and tt- junctions 8,9 . In a pio- 
neering workifi Bulaevskii et al. demonstrated that the 
spontaneous Josephson vortex (JV) carrying flux <5>o/2 
appears at the boundary between 0- and n- JJs. The 
structure of such semifiuxon has been studied in more 
detail in Refs. HlTl 

In the present work we study the properties of the pe- 
riodic array of 0- and tt- JJs. We show that depending 
on the ratio of lengths of 0- and tt- JJs, such an array can 
be either in the homogeneous or modulated state. The 
second order phase transition between these states takes 
place when the length mismatch between 0- and tt- JJs 
is small. At the transition point the screening length of 
the magnetic field diverges. In the modulated states the 
average phase difference </?o can take any value between 
— 7r and 7r. Further, we call such systems tp-junctions. 
Such structures were first predicted by Mints^ in the 
case of randomly alternating 0- and tt- JJs. We study 
the properties of (/?-j unctions: effective Josephson cur- 
rent, Josephson length, and Fraunhofer- like oscillations 
of the critical current. Magnetic properties of long ex- 
junctions are determined by two types of Josephson vor- 
tices carrying partial fluxes $o</?o/ 7r an d $o( 7r — Po)/ 7 ? 
(see alsoiiii) . We find analytically the shapes of these 
unusual Josephson vortices. We also demonstrate that 
at the boundary between (p- and usual JJs a Josephson 
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FIG. 1: Phase diagram of the periodic alternating 0-tt junc- 
tion structure. Inset shows dependence of the equilibrium 
phase on the length mismatch for dorfvr = Aj 

vortex appears carrying a partial flux ^qipq/2tt. 

Consider a periodic structure composed of alternating 
0- and 7r-JJs of lengths do and d v respectively. The en- 
ergy per period of such a structure is given by 

'-s/l($(S)' +i -*)* 

where j c is the Josephson current density, Xj is the 
Josephson length of the individual junctions, X 2 j = 
c<f>o/(87r 2 tj c ), and t = to+Ai + A2 is the effective junction 
thickness^. For simplicity, we focus on the case when the 
critical current densities j c are the same in 0- and 7r-JJs. 
Generalization for different critical currents is straight- 
forward. Possible realization of such system is the zigzag 
junction between high-T c and conventional superconduc- 
tors as well as S-F-S junctions with periodically modu- 
lated thickness of the ferromagnetic interlayer. In the 
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case of the zigzag structure, x is the coordinate along the 
zigzag boundary. The ground state phase distribution is 
determined by the equation 



(2) 



A , 2 in 7r-JJs. At 



where j(x) = — A f 2 in 0-JJs and j(x) 
the boundaries, the phase and its derivative have to be 
continuous. Periodicity of </>(x) implies that the solution 
is symmetric with respect to x = —do/2 and d n /2, i.e., 
d(j)/ dx\ x= _ do f 2 — d4>/dx\ x —d^/2 = 0. Eq. always has 
homogeneous solutions = and = 7r. However these 
solutions give the ground state only in some range of 
the ratio do/d^. In general, Eq. also allows for the 
inhomogeneous solution. Consider the case d^ < do an d 
weak modulation around the = state. In this case 
Eq. J21 can be linearized and its solution is given by 
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—do < x < 0, 

< x < dir. 



(3) 



Matching the logarithmic derivative dhiip/dx of these 
solutions at x — 0, we obtain the following condition 



tanh(<io/2Aj) = tan(d w /2Aj), 



(4) 



for the onset of the modulated solution. The energy 
analysis shows that the uniform — solution is favor- 
able at tanh (d /2Aj) > tan (d v /2Xj). In the opposite 
case, dir > do, the uniform — tt solution is favorable 
at tanh (c? T /2Aj) > tan (do/2Xj). In the case of a fi- 
nite junction split into 0- and tt- pieces, condition Q 
was first derived in Ref. ITol These results are summa- 
rized in the phase diagram shown in Fig. ^ In the limit 
do 3* Aj the condition for the modulated solution is given 
by dir > (w/2)Xj. In the case do,^ <C Aj the region of 
the modulated solution is given by 
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i.e., the modulated solution exists only for a very small 
length mismatch. Since Xj is temperature dependent, 
the transition into the modulated state may occur with 
decreasing temperature. 

In the case do = d^ = d, the phase distribution is al- 
ways modulated. Lets find this distribution. From sym- 
metry 4>{—x) = tt — 4>{x), in particular, (f)(0) = tt/2. The 
first integral of Eq. J21 ) for x < is given by 



A 2 , 



dx 



cos < 
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with <po = 4>(—d/2). The solution is given by 



a/cos 0o — COS 0' ' 



(6) 
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FIG. 2: Le/f panel: dependence of the minimum phase on the 
length of junction for do = d^. Right panel: shape of phase 
modulation at different junctions sizes 



where <po is determined by the condition x(4>o) = —d/2. 
Dependence 4>o(d/\j) and shapes of the phase variation 
4>{x) for different values of d/Xj are presented in Fig. [21 
For small length d <C Xj, it is easy to find that the phase 
distribution for — d < x < is 4>(x) = tt/2 + x(x + 
d)l (8Aj), i.e., it varies weakly around 7r/2. 

Consider a structure composed of short alternating 
and tt- JJs with a small length mismatch, \d$ — d n \ -C do- 
In this case the function j (x) in Eq. @ oscillates rapidly 
and we can present <j>(x) as <f>(x) — <p(x)+£(x), where the 
function ip(x) varies slowly over the distance do and 

0. In the limit do j7r < Xj, 



oscillates rapidly with £(x) 
the oscillating term is small (as one can see from the 
right panel of Fig. |21 even at do,d w = Xj, < 
0.13 <C 1 and (ii) the oscillating part is smaller than 
the average phase). In this case one can use a coarse- 
grained approximation (see, e.g., Ref. meaning that 
we can average out the rapidly oscillating function £(x) 
and derive the equation for the slowly varying function 
y>(x) only: 
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(do +dir)X 2 j 



sin (</>), (8) 



Here the first term in the right hand side appears due to 
current modulations and the second term is proportional 
to the average Josephson current. A similar equation 
was derived in Ref. 13 . The neglected terms in Eq. © are 
smaller by a factor of (do/Xj) 2 . Without the external 
magnetic field, when condition 101 is fulfilled, the equi- 
librium phase difference ipo is: 



cos (fio = I 
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dir) 



d Q dir (d + d^) 



(9) 



and for |/| < 1, the average phase difference lies in the 
region < \(po\ < 7T,i- e. we have the realization of a J J 
with arbitrary ground state phase difference ( (p- junction). 
For |/| > 1 homogeneous 0- or tt- phase is realized. An 
example of the dependence of the equilibrium phase tpo 
on length mismatch is shown in the inset of Fig. ^ 
The energy in terms of the coarse-grained phase is 
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minimum value of A c ff, A ff = 2d\j yfiido — d^) is 
reached at / = 2 

Exactly at the transition point, / = 1, screening is 
nonlinear and the phase variation obeys equation 



2X 2 v d 2 tp/dx 2 = p 3 



(15) 



The general solution of this equation is given by tp — 
±2X v /(x + C). Using the phase-field relation H = 
($o/2irt)d(p/dx, we derive from this solution that the 
magnetic field decays inside the structure as 



FIG. 3: Left panel: current-phase relationships for different 
values of the equilibrium phase tpo (see Eq. (|1 ip 'I . Curves 
are marked by the value of tpo- Unit of current density j in 
this plot is d ! j c /(l2X 2 J ). Right panel: dependencies of critical 
current of small 0, tt- structure on the magnetic flux $ through 
the structure, for different values of tpo- The current unit at 
the vertical axis is jj Sdo/(12Xj), where S is the cross section 
area of the structure. 



given by 
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1 ( dp\ 2 (cos ip — f) 2 



dx 



2X% 



dx, (10) 



with A 2 = \2\^jd\ 3> Aj. The coarse-grained Josephson 
current density flowing through the array is given by 

4i 
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sin {tp) (/ - cos ip) . 
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This current-phase relation is quite peculiar (see left 
panel in Fig. EJ). The current has two maxima 
and two minima, which are achieved at cos^i = 



/ ± y/f 2 + 8 J /4. The absolute maximum 
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is achieved at tp — p-. 

Consider screening of the weak magnetic field by the 
periodic structure. In the modulated phase linearizing 
Eq. JSJl near tpo, tp = po + tp, we obtain a linear equation 
for the phase variation tp 



d tp sin po _ 
dx 2 



-tp, 



(12) 



One can see that the effective screening length 



A cff = A^/sin^ = X v /y/l - p, for |/| < 1 (13) 

diverges at the transition point. In the homogeneous 
phase, |/| > 1, similar derivation gives 



Aeff = A^/Vl/l-1, for |/| > 1. 



(14) 



Because X V (T) diverges at T — > T c faster than /(T), 
A ff (T) is nonmonotonic above the transition point. The 



H(x) = H 1 + x/ 



Ay$Q 
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where Hq is the external magnetic field. It is interesting 
to note that the field tail at large x does not depend on 
the applied field 
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for x 3> 
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In addition to the usual 27r-degeneracy, in the case / < 
1, the ground state energy is also degenerate with respect 
to sign change of the equilibrium phase tpo — > — po- This 
degeneracy leads to the appearance of two new kinds of 
solitons in which the phase sweeps either from —po to 
tpo or from po to 2n — pq. Rewriting Eq. JHJl in a more 
convenient form 
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sin (p) (cos tp — cos tp) , 



we derive the shapes of the two solitons 
pi(x) — 2arctan ^tanh ^ — — ^° ^ tan ("7^")^ ; (18 



a 



tp 2 (x) — tt + 2 arctan 



/ tanh(xsin p /2X v ) 



\ tan (tpo/2) 
The soliton energies per unit length are given by 



(18b) 



£i = e./o(Aj/4A0) (Isin^ol - tpocostpo) , (19a) 
£2 = e./o(Aj/4A0) (|sin^ | + (tt - tp )cos<p ) , (19b) 

where e,jo = 4$oj c Aj/(7rc) is the energy of a single soli- 
ton in a uniform J J. At the transition point, / = 1, the 
first soliton vanishes and the second one acquires the fol- 
lowing shape, p 2 (x) = tt + 2 arctan (x/X v ). Even though 
for do > d v the first energy is smaller than the second 
one, both solitons play a role in the magnetic properties 
of the (/^-junction. This is because without the second 
solitons the system can not contain two or more first 
solitons of the same sign due to topological constrains. 
In the external magnetic field H , the energy of the sys- 
tem with a small number of solitons of both kinds, N\ 
and N%, can be written as 

E = s 1 N 1+ s 2 N 2 -(n 1 ^+N 2 ^-)^, (20) 

\ TT TT I 47T 
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and has to be supplemented by the topological constrain 
\Ni — N 2 \ < 1. The penetration scenario is determined 
by the ratio of the soliton energy e a to its flux $ Q . Simple 
analysis shows that for do > d v these ratios are arranged 
in the order £i/$i < (si +£2)/ < i > o < £2/^2- These in- 
equalities imply that a single first soliton appears at the 
field Hi = 47T 2 £i/</?o < 5 ) o- The soliton lattice, composed 
of alternating solitons of two types, penetrates into the 
structure at higher field H c \ = 47r(ei + 62) /§o- In the 
case do — d v , the fields Hi and H c i coincide. 

An interesting possibility to have an equilibrium ip 
Josephson vortex is realized at the boundary between an 
alternating structure of 0- and n- JJs and usual 0-JJ. We 
assume that the alternating structure occupies the region 
x < 0, while the 0-JJ with the Josephson penetration 
length A occupies a positive semi-axis x > 0. In this case 
the phase difference must vary continuously from ipo for 
x — ► — 00 to zero for x — > +00. Due to the condition 
X v ^S> A ~ Aj, the value tp(0) -€i 1, which means that the 
boundary vortex is almost completely localized in region 
x < and it is equivalent to half of the first soliton (|18a|l . 
The metastable n — 2ipo boundary vortex is also possible 
and in the lim it X v 3> A its shape corresponds to the 
second soliton (|18bjl . 

Consider a structure with the total length L smaller 
that the screening length A e ff in an external magnetic 
field H, smaller than <&o/(do + d^t. In this case the 
coarse-grained approximation for the phase is justified 
and we can use the current-phase relation where 
the coarse-grained phase depends on the coordinate as 
ip(x) — hx + (3 with h — 2ntH / $0 and /3 is an arbitrary 
phase shift. The total current per unit thickness flowing 
through the structure is given by 

T(R\ ^ Tw (t ■ fh sin(2(fe a ; + /3)) 
/ (P) = Y2A J / 1 f sm ( hx +P) 2 

J -L/2 

Calculating the integral, we reduce it to the form 

lift) — — sm (fi) sm V (/ — cos W) cos v) 7 (21) 
V 



where I Q = j c Ldl/(12X 2 j) and r) = hL/2 = 7r$/$ and 
<f> = tLH is the total flux through the structure. The 
current reaches a maximum at 



cos 09) = C0S ^~^ C ° S2 ^ +8C0S2(?7) . (22) 
v ' 4 cos (77) v ' 



Shapes of the field dependencies of the critical current 
at different values of ipo are plotted in the right panel of 
Fig. 13 These dependencies differ significantly from the 
Fraunhofer dependence in usual JJs. They have large a 
component with half of the main period. At ipo = ir/2 
the dependence has the Fraunhofer shape but the pe- 
riod is two times smaller than in usual case. Note that 
the observed in the alternating 0-ir junction structure 8 
finite critical current at zero magnetic field is naturally 
obtained in the framework of our analysis. The theoret- 
ical analysis in Ref. 0, giving zero current at H = 0, is 
incomplete, because it does not take into account the cur- 
rent term oc sm(2ip) coming from the rapidly oscillating 
phase. 

In conclusion, the possibility to realize the transition 
into the ip— junction state with decreasing the tempera- 
ture from T c , may be very helpful for experimental ver- 
ification of the predicted effects. In particular, the ob- 
servation of the striking non-monotonous variation of the 
screening length with temperature would provide an un- 
ambiguous proof of such a transition. Also, the studies 
of the fine structure of the critical current dependence vs 
magnetic field and the observation of a periodicity two 
times smaller than that expected for standard J J, could 
be of considerable interest. Finally, note that the SQUID 
microscope may directly probe the specific shapes of the 
partial flux vortices. We would like to thank A. Lopatin 
and W. K. Kwok for critical reading of the manuscript. 
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